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We discuss the problems that arise when one wishes to extend the existing general methods of 
computing radiative corrections to QCD jet cross sections to beyond next-to-leading order. Then 
we present a subtraction scheme that can be defined at any order in perturbation theory. We give 
a status report of the implementation of the method to computing jet cross sections in electron- 
positron annihilation at the next-to-next-to-leading order accuracy. 
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Accurate predictions of QCD jet cross sections require the computation of radiative corrections 
at least at next-to-leading order (NLO) accuracy, but in some cases also at higher order. The physi- 
cal cases when computations at the next-to-next-to-leading order (NNLO) are important have been 
discussed extensively in the literature [JI]]. There are also some less standard motivations. Although 
several general methods exist for computations at NLO, these become very time-consuming with 
increasing number of final-state partons. Many talks at this conference showed impressive progress 
in computing one-loop corrections to amplitudes of multileg processes [0], therefore, we should 
consider the question wether we can compute the real radiation corrections fast enough. Further- 
more, the thorough understanding of computing NNLO corrections may help the combination of 
parton showers and NLO calculations. 

The perturbative expansion of any jet cross section can formally be written as a = a LO + 
a NLO _|_ ^nnlo _|_ Let us consider e + e~ —> m jet production, when a LO is the integral of the 
fully exclusive Born cross section over the available phase space defined by the jet function J m , 

a LO = [ dalJ m = [ d$ m \^ m 0) \ 2 J m . (1) 



There are two contributions to the NLO correction. We have to consider the fully exclusive cross 
section da R for producing m + 1 partons and the one-loop correction da v to the production of m 
partons, 

a NLO = / da R 7 m+ i + / da v 7 m = /d0 m+1 \Jt [ *\ x \ 2 J m+l + /d0 m 2Re(^ 1) \JK^)J m . (2) 

Jm+\ Jm J J 

These two contributions are separately divergent in d = 4 dimensions although their sum is finite 
for infrared safe observables. We assume that ultraviolet renormalization has been carried out, 
so the divergences are purely of infrared origin and are regularized by defining the integrals in 
d = 4 — 2e dimensions. 

There are several general methods of computing the finite NLO correction. Most of these rely 
on the same principles, namely one defines approximate cross section da A which regularizes the 
real correction in d dimensions in all its infrared singular limits that lead to poles in e, so the cross 
section 

<£$=[ [(da R ) e=0 / m+ i-(da A ) £=0 7 m ] (3) 

Jm+\ 

is finite. 1 The approximate cross section is constructed using the universal soft- and collinear 
factorization properties of QCD matrix elements (we use the colour-state notation [J3j] and also 
some notation introduced in Ref. [Q]), 2 

S r \^i%(p r ,...)\ 2 oc^^(^l \...)\TrT k \^ m °\...)), (4) 

; ^ SirSkr 

i+k 



'The subtraction terms in these equations are symbolic in the sense that these are actually sums of different terms. 
The jet function depends on different momenta in each of these terms, the exact set of momenta for each term can be 
found in Ref. JI|]. 

We drop some numerical factors in order to keep the expressions as simple as possible, as only the structure of 
these formulae is relevant for the discussion. 
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d(°)i sA°), 



C ir \^^(p i , Pr ,...)\ 2 o,-(^r(pir,...)\Pr\^rn > (Pir,---))- 



Si, 



(5) 



These factorization formulae allow for such a construction that the integration over the phase space 
of the unresolved partem can be computed independently of the jet function, leading to 



^da A = da,>/(£), 
where 1(e) is an operator in the colour space with universal pole part, 



(6) 



a s 



'<">-£E 



1 1 Vr t (4 ¥ 2 



+ 0(£ C 



(7) 



da v +/da A 



(8) 



e=0 



This pole part is equal, but opposite in sign to the pole part of the virtual correction, so that the 
ra-parton integral 

Jm I J\ 

is finite. Therefore, the sum of the two finite contributions a,^ LO and G^ffi is equal to a 1 ^ . 

In constructing the approximate cross section special care is needed to avoid double subteac- 
tions in the regions where a soft parton becomes also collinear to another hard parton. At the NLO 
accuracy, the overlap of the soft- and collinear limits can easily be identified to be the collinear 
limit of the soft factorization formula [|J]. However, disentangling the multiple unresolved limits 
at higher orders, when multiple soft-, collinear- and soft-collinear limits overlap in a complicated 
way, is far more cumbersome [Q], This calls for a simple and systematic procedure. 

In a physical gauge the collinear singularities are due to the collinear splitting of an exter- 
nal parton [Q], The overall colour structure of the event does not change, the splitting is entirely 
described by the Altarelli-Parisi functions which are a product of colour factors and a kinemati- 
cal function describing the collinear kinematics of the splitting. Thus, if we want to identify the 
collinear contributions in the soft factorization formulae to any order in perturbation theory, we 
can use the following simple procedure: (i) employ the soft insertion rules [H [7|] to obtain the usual 
expression 



,M„v 



in in o r* v 

S^l l {p n ...)\^Y J Y J Y j e,{p r )e;{ Pr )-^^{^\...)\T v T k 



i=U=lhel. 



Sir^kr 



(0) 



with 



s ir = 2p rPr and £ s,( Pr )s;( Pr ) = -g^ v + Pr - +P> ^ 

hel. P r ' n 



(..-)), (9) 



(10) 



(ii) choose Coulomb gauge (n* 1 = Q^ 1 — Pr Q 2 /s r Q, SrQ = ZprPo) to identify the collinear contri- 
bution in the colour diagonal terms 



S r K U i(^,...)| 2 oc£ 



1=1 




(^i 0) (...)\T r T k \^i 0) (...)) 



(11) 
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(iii) define momentum fractions in the Sudakov parametrization of momenta pf and pr being 



collinear as zi 



>iQ 



SiQ+S,-Q 



, so that the colour-diagonal terms become equal to the collinear limit of the 



soft factorization formula. Then the pure soft contributions are given by 



sn^SiCpr,-)! 2 ' 



I 

1=1 



i 



sik 2s{q 2skQ 
: l^ I 



(^n\...) 



■T k \^i 0) (...)) 



(12) 

We checked explicitly that this procedure leads to non-overlapping factorization formulae that de- 
scribe the analytic behaviour of the squared matrix elements in any IR limit at the NNLO accuracy 
[B]. Furthermore, the factorization formula in the purely soft limit is independent of the helicity of 
the soft gluon. This allows for the definition of approximate cross sections for real radiation with 
fixed helicities, and thus for Monte Carlo summation over the helicities in NLO computations. The 
Monte Carlo summation over the helicities in computations at LO was found very useful for saving 
CPU time[|]. 

The physical motivation for higher accuracy and the success of the subtraction schemes at 
NLO lead one to consider the extension of the subtraction method to NNLO, when three terms 
contribute: the double-real, the real-virtual and the double-virtual cross sections, 



,-NNLO _ ,_RR , _RV , _VV _ / A„RR j , / Ht RV J _l_ / Ht VV J 

Im 



J m+2 J m-\ 



_ rl j m +i+ I da,; v 7 m . (13) 

lm+2 Jm+\ Jm 

The necessary ingredients for constructing approximate cross sections, namely (i) the tree level 
three-parton splitting functions [10] and double soft gg and qq currents [{7|, [ll|] and (ii) the one-loop 
two-parton splitting functions [ |l2| ] and soft-gluon current [13] that is, the infrared (IR) structures 
of the three contributions at NNLO have been known for some time. The difficulty of using the 
multiple infrared factorization formulae for cunstructing the approximate cross sections is amply 
demonstrated by the slow progress in setting up a subtraction scheme. Other approaches to com- 
puting NNLO corrections have been more successful. 

The antennae subtraction method [|14]] uses complete squared matrix elements instead of the 
IR structure and the first complete computation of NNLO corrections to three-jet production in 



electron-positron annihilation has been reported to this conference [15]. For processes involving 



massive particles and/or simple kinematics, direct numerical evaluation of the coefficients in the 
Laurent expansion of the three contributions (based on sector decomposition) lead to the complete 



NNLO corrections of Higgs- [ |16| ] and vector-boson production in hadron collisions Q17Q 
The reorganization of the NNLO contributions into three finite cross sections, 



.NNLO _ ^NNLO , ^NNLO , ^NNLO 



'tn+2 



+ °m+l + O, 



m i 



(14) 



is governed by the jet function as follows: 



and 



-W2 



r NNLO 



Jm+2 <- 



r NNLO _ / J j/T-RR J H,r RR > A 2 J 



J m+2 J m 



'm+\ 



m+1 



{(dcCi+/^* Xi !.4,, 



, RR,A! T ,„RR,Ap 7 \\ 

dCJ m+2 J m+l-aO m+2 Jm) j , 



*<l? l + {jteT^ A > 



(15) 
(16) 



r NNLO 



daZ + / (da m+2 2 -^ m+2 



do»" ) + 



da-f 1+ (y id c 2 AiVv 



(17) 
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Let us concentrate on Eq. (|16|). The construction of da m+ \ ' that regularizes the kinemati- 
cal singularities of do^Yj in the singly-unresolved regions is straightforward, but the difference 



da m+1 7 m+ i-da m+1 7, 



is infinite. To make it finite, we have to subtract the universal pole 
part, given by Eqs. (|6])-(|8|), too. The latter however, does not obey universal collinear factorization. 
Due to coherent soft-gluon emission from unresolved partons only the sum {-^ m l\\{T j-Tk + T r - 
Tk)\^ml\) factorizes in the collinear limit (Tj r = T j + T r ), 

Cjri^AiTj-Tt + Tr-T^^^oc—^lTjr-TtP^^). (18) 

Sj r 

This factorization is violated by the factors s~[ k e /e 2 at 0(e), which was also noticed in Ref. fll9|], 
where it was shown that the terms that violate factorization are known to give vanishing con- 
tribution after integration. However, if one insists on defining fully local subtractions, which is 
important for numerical stability and reducing CPU time, then the use of properly defined new ap- 
proximate cross sections is necessary. The complete subtraction scheme at NNLO, based on these 



new approximate cross sections is defined in Refs. [18, ^Op.We employed this subtraction scheme 
for computing the finite cross sections o^S? and d^ 1 ^ of the C-parameter and thrust distribu- 
tions in electron-positron annihilation. In order to have the complete physical prediction we also 
have to compute a^^ , which requires the integration of the subtraction terms over the singly- 
and doubly-unresolved factorized phase spaces. We used standard techniques fractioning, sector 



decomposition [ |21[ ] and residuum subtraction to find the Laurent expansion of the one-particle 
integrals in 

"'da^ + ( /da* R ^ A 'l. (19) 



*Cf 2 -d<Cf 2 )- (2°) 



We expect that the same techniques can also be employed for the computation of the coefficients 
in the £-expansion of the two-particle integral 

This work is in progress. 
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